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Starting from an effective action for the order parameter field, we derive a 
coupled set of generalized hydrodynamic equations for a Bose condensate in 
an optical lattice at finite temperatures. Using the linearized hydrodynamic 
equations, we study the microscopic mechanism of the Landau instability due 
to the collisional damping process between condensate and noncondensate 
atoms. It is shown that the Landau criterion of the superfluidity for the 
uniform system is modified due to the presence of the periodic optical lattice 
potential. 

1. Introduction 

Recently Bose condensed gases in an optical lattice are extensively stud- 
ied both experimentally and theoretically. Much of attention is attracted 
to the breakdown of the superfluidity in a condensate. A lot of works 
have discussed the stability condition of the superfluidity by using the zero- 
temperature Gross-Pitaevskii equation. 1-4 As emphasizing in our previous 
works, 5 it is crucial to explicitly include the degree of freedom of noncon- 
densate atoms into the theory to understand the Landau instability from 
the microscopic point of view. 

In the present paper, by means of an effective action and an associated a 
coupled set of generalized hydrodynamic equations, we study the microscopic 
mechanism of the breakdown of the superfluidity having experiments such 
as reported in Ref. 6 in mind. 
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2. Generalized Hydrodynamic Equations in an Optical Lattice 

We consider the Bose condensate at finite temperatures trapped in a 
one-dimensional optical lattice in addition to the harmonic potential 



Kxt(r) = —{oo 2 x x 2 + uj 2 y y 2 + lo 2 z 2 ) + sE R cos 2 ( -z ) , (1) 



where uj x , uj v , u z are the frequencies of the harmonic trap potential, s is 
the strength of the lattice potential, En is the recoil energy, and d is a 
lattice constant. In the external potential represented by Eq. (1), the Bose 
condensates are described by the following action 



S = dv dt 6*(rA 



2gh{r, t) 



+ Y J drdv' J dtdt'(t)*{r,t)F{r,t;r',t')(l){v',t'), (2) 

where 4>(r, t) is a complex field that represents the degree of freedom of 
condensate atoms, n(r,t) = (^(r, t)ip(r, t)} is the noncondensate density 
with tp(r,t) being the noncondensate field operator, and F(r, t;r',t') = 
G(r, t; r' , t')G(r, t; r', t')G(r',t'; r, t) is defined in terms of nonequilibrium Green's 
functions of the noncondensate atoms G(r, t;r',t') = — i(ip(r, t)ip^(r', t')). In 
Eq. (2), g is a coupling constant which is related to the s-wave scattering 
length of atoms. We note that the index c of integral in Eq. (2) means 
that time integrations should be performed on Keldysh contour path. 7 The 
derivation and detailed discussion of the action (2) will be presented else- 
where. 8 

Considering a relatively deep lattice potential with large s, we introduce 
the ansatz (tight-bindeing approximation) 

<f>(r,t) = Y^Mz)Mr±,t)e i6l{Tx ' t) , (3) 
i 

for the order parameter field, where 6i(r±,t) is the phase of the order pa- 
rameter at the Zth lattice cite, and wi and <pi are real functions. Here we 
assume the periodicity for wi{z) = wq{z — Id), where wo is localized near 
z = 0. Substituting the ansatz (3) into the action (2) and performing a 
smoothing procedure from discrete variables to continuous ones following 
Ref.9, we obtain after some algebra 



S = dt dr 



yjn c (r,t)ih^-^jn c (r,t) - %n c (r,t) d 
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m(r,t) 



+ J dt drn c (r, t) cos d 



Oz 



^JdtJ dr^/n c (r,t)|vi^/n c (r,t) - yf^rj) (V ± 9(r,t)Y 



2V ±v /n c (r,t) • V±0(r,t) + yjn c (r, t)V 2 ± 0(r, t) 
Jdtjdr V ho (r)n c (r,t) - | J dt J dr n 2 c (r,t) 

ae(r,t) 



J J dt J dr n c (r, t) cos yl 
2g J dt J dr n c (r,t)h(r,t) 



9 2 



dz 



dtdt' I drdrJn c {r,t)e- i9{r ' t) F(r,t;r',t')Jn c (r',t')e 



(r,t) 



(4) 



where g = gd J dzw^z) is the effective renormalized coupling constant, 
Vho(r) is the harmonic potential in Eq. (1), and n c (r,t) = cpf(r±,t)/d and 
9(r,t) = Oi(r±,t) represent the density and phase after the smoothing pro- 
cedure. In Eq. (4), we defined the following quantities: 



J 
J 



J dz wi(z) 



tf_SP_ 
2m dz 2 



V op t{z) 



Agd / dzwf(z)wi + i(z) 



(5) 
(6) 



where V op t(z) is the optical lattice potential in Eq. (1). 

The hydrodynamic equations are derived by imposing the stationary 
condition on the smoothed effective action (4) with respect to arbitrary 
variations of the density n c (r,t) and of the phase 6(r,t). After performing 
the Markovian approximation and the gradient expansion for the non-local 
term related to F(r, t; r', t'), we obtain 



m — h V 

at 



~2^ Vz{r,t) + 2" V ^ r '^ + ho(r) 



+ gn c (r,t) + 2gh(r,t) 
dn c (r,t) d 



0, 



(7) 



dt dz 



rn 



m. 



n c (r,t)v z (r,t) 



+ V± (n c (r, i)v ± (r, t)) = -T 12 (r, t). 

(8) 
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where Ti2(r,t) = 2n c (r, t)R(r, t) with 



* dki dk.2 dk.3 



2h J (2vr)3 (2tt)3 (2tt) 3 
x 5(uj c +&1-U2- a3 3 )<5(k c + ki - k 2 - k 3 ) 

x I [1 + /(ki, r, t)] [1 + /(k 2 , r, i)]/(k 3j r, i) 
-/(k 1 ,r,t)/(k 2 ,r,t)[l + /(k 3> r,t)]} ) 



(9) 



where p c = ftk c and e c = hco c is a condensate momentum and energy, / is 
the nonequilibrium distribution function of the noncondensate, and Tiooi = 
h 2 k 2 /2m + 2g(n c (r,t) +n(r,t)) is a noncondensate energy. The condensate 
chemical potential and energy are defined by 



Mc(r,*) 
e c (r,*) 



gn c (r,t) + 2gn(r,t), 



(10) 



A*c(r,t) + y^(r,t) + y^M) + y f — J t£(r,t). (11) 



In deriving above hydrodynamic equations (7) and (8), we have expanded the 

gradient term of the phase regarding the term as small, which is the case in 

the study of small amplitude oscillations, and introduced the effective mass 

-1 



and the chemical potential effective mass, m* = h 2 /d 2 J — Jn c (r, t) 



and 



m* = h 2 /d 2 J — 2Jn c (r,t) , respectively. We neglected the n c depen- 
dence of m*, quantum pressure term, and VR(r,t), in deriving Eqs. (7) 
and (8). The coupled set of hydrodynamic equations (7) and (8) are useful 
starting point for discussing the collective oscillations of the condensate in 
an optical lattice at finite temperatures. 10 > n > 6 



3. Instability due to the Collisional Damping Process 

In this section, we demonstrate the useful application of the hydrody- 
namic equations (7) and (8) derived in Sec. 2. by investigating the stability 
of the superfluidity. We consider the system in which the condensates are 
moving with a constant velocity 6 and neglect the confining trapping poten- 
tial. 

Nonequilibrium distribution functions in Eq. (9) should be determined 
by solving a kind of kinetic equation. Instead of solving a kinetic equation, 
we approximate the nonequilibrium distribution function by the static Bose 
distribution function /(p,r) = e p^ p -n)_ v 12 where e p = Hlu p = ^ Fi +2g{n c + 
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n) and jl is the chemical potential of the noncondensate atoms. Expanding 
e c (r, t) around equilibrium, one can linearize r^r, t) as 

-g 2 f dki dk.2 dk.3 



,r 12 ,M) = f/ 



r<5(w Cj0 + Wl ~U 2 -^3) 



(2tt) 3 (2tt) 3 (2vr) 3 
x 5(k c n + k x - k 2 - k 3 )[l + /(k l5 r)]/(k 2 , r)/(k 3 , r) 



x p 
1 



g5n c + mvj_ • £vj_ + m 



m 
rri' 



v z ,o$v z 



/' / 



r'(r,t) 



gbn c + mv i)0 • 5v± + m — f 2j0 ^ 



(12) 



where we define the relaxation time r'(r, t) by 
1 7T<? 2 /? /" dki dk.2 dk.3 



r>(r,t) 



2T-I 



(2tt) 3 (2^) 3 (2vr) 3 
x 5(a; Ci o + ujx - uj 2 - uj 3 )5(k cfi + ki - k 2 - k 3 ) 

Finally we obtain the linearized hydrodynamic equations 



(13) 



^M + V,.Ko(^v c (r,t)] + | 

+ V ± -[5n c (r,t)v c>0 (r)] + — 

oz 



111 



rrr 



—5n c (r,t)v zfi (r) 



-sr 12 (r,t), 



(14) 



9 r , \ _ 

m—dv c (r,t) = -V 



^v_L )0 (r) • 5vj_(r,t) + m — v z ,o(r)^ 2 (r,t) 



m; 



+ #<5n c (r,f) 



(15) 



Looking for a plane-wave solution such as exp[i(k • r — tut)] with k = (0, 0, k), 
we obtain uj = v^ • k + c*k = £1 in the limit 1/r' — > 0, where c* = 



gn c fi/m(m/m*) is the sound velocity modified due to the lattice potential 
and v^ = (m/mjjvcfi is the chemical potential group velocity. 3,4 When we 
include 1/r', the dispersion relation is modified as u = $l — iT. To first order 
in 1/r', we obtain the damping term 



2r' \ m* c* , 



(16) 



where v Cj o and k are normal vectors for v Cj o and k, respectively. From the 
expression (16) for T, it can be seen that the instability occurs for the mode 
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with v Cj o • k = — 1 when > c*, 2-4 since the damping term T changes its 
sign. 



4. Summary 

In summary, we have investigated the stability of the superfluidity of 
the condensates in an optical lattice by using the generalized hydrodynamic 
equations. We derived the Landau criterion microscopically and showed 
that the criterion is modified due to the lattice potential as expected. In 
the present paper, while we only consider the collisional damping process 
as a destabilizing process, one can also derive the Landau instability due 
to the Landau damping process, which is also important in the collision- 
less regime of interest, from the generalized hydrodynamic equations. 8 The 
generalized hydrodynamic equations derived in the present paper can be 
used for analyzing the dipole oscillations and its damping in the presence of 
the noncondensate atoms such as observed in the experiment by Florence 
group. 10 
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